(19) 



3 



Europaisches Patentamt 
European Patent Office 
Office europeen des brevets 



(12) 



(11) EP 0 818 131 B1 

EUROPEAN PATENT SPECIFICATION 



(45) Date of publication and mention 
of the grant of the patent: 
13.10.2004 Bulletin 2004/42 

(21) Application number: 96907215.6 

(22) Date of filing: 28.03.1996 



(51) Intel J: H05G 1/02, G21K 1/02, 
G03B 42/02, A61B 6/00, 
G01N 23/04 

(86) International application number: 
PCT/AU 1996/0001 78 

(87) International publication number: 

WO 1996/031098 (03.10.1996 Gazette 1996/44) 



(54) SIMPLIFIED CONDITIONS AND CONFIGURATIONS FOR PHASE-CONTRAST IMAGING WITH 
HARD X-RAYS 

VEREINFACHTE BEDINGUNGEN UND KONFIGURATIONEN FUR 
PHASENKONTRAST-BILDERZEUGUNG MIT HARTEN RONTGENSTRAHLEN 

CONDITIONS ET CONFIGURATIONS SIMPLIFIES POUR IMAGERIE A CONTRASTE DE PHASE 
A RAYONS X DURS 



(84) Designated Contracting States: 

AT BE CH DE DK ES Fl FR GB GR IE IT LI LU MC 
NL PT SE 

(30) Priority: 28.03.1995 AU PN201295 

(43) Date of publication of application: 
14.01.1998 Bulletin 1998/03 

(73) Proprietor: X-Ray Technologies Pty Ltd 
Melbourne, Victoria 3000 (AU) 

(72) Inventor: WILKINS, Stephen, William 
Blackburn, VIC 3130 (AU) 

(74) Representative: Charlton, Peter John 
Elkington and Fife LLP, 
Prospect House 

8 Pembroke Road 

Sevenoaks, Kent TN13 1XR (GB) 



CD 

CO 

CO 
00 

Q_ 
LU 



(56) References cited: 
EP-A- 0 274 155 
EP-A- 0 599 361 
US-A- 4 200 354 
US-A-5 016 267 



EP-A- 0 466 047 
WO-A-95/05725 
US-A- 4 953 188 
US-A- 5 319 694 



T.J. DAVIS ET AL.: "PHASE-CONTRAST 
IMAGING OF WEAKLY ABSORBING 
MATERIALS USING HARD X-RAYS M NATURE., 
vol. 373, 16 February 1995, LONDON GB, pages 
595-598, XP002033451 

JOURNAL OF MICROSCOPY, Vol. 161, Pt. 3, 

March 1991, SHINOHARA et al., "Radiation 

Damage in Soft X-Ray Microscopy of Live 

Mammalian Cells", pp. 463-72. 

OPTICS (Hecht/Zajac), 1973, pp. 449, 450 and 

474-479. 

WHITE V. AND CERRINA F.: 'METAL-LESS 
X-RAY SPHASE-SHIFT MASKS FOR 
NANOLITHOGRAPHY' JOURNAL OF VACUUM 
SCIENCE & TECHNOLOGY vol. B10, no. 6, 
November 1992 - December 1992, NEW YORK, 
US, pages 3141 -3144 



Note: Within nine months from the publication of the mention of the grant of the European patent, any person may give 
notice to the European Patent Office of opposition to the European patent granted. Notice of opposition shall be filed in 
a written reasoned statement. It shall not be deemed to have been filed until the opposition fee has been paid. (Art. 
99(1 ) European Patent Convention). 



Printed by Jouve, 75001 PARIS (FR) 



BEST AVAILABLE COPY 



EP0 818 131 B1 

Description 

[0001] This invention relates generally to the observation of a structural feature of an object utilising penetrating 
radiation such as x-rays. More particularly, the invention relates to x-ray contrast recordal and imaging of internal 
5 boundary features. 

[0002] The present applicant's international patent publication WO95/05725 (PCT/AU94/00480) discloses various 
configurations and conditions suitable for differential phase-contrast imaging using hard x-rays. Other disclosures are 
to be found in Soviet patent 1402871 and US patent 5319694. It is desired that relatively simpler conditions and con- 
figurations more closely related, at least in some embodiments, to traditional methods of absorption-contrast radiog- 
10 raphy, may be utilised for differential phase-contrast imaging with hard x-rays. 

[0003] The present invention provides a method of deriving a phase-contrast record of a boundary as set out in 
claims 1 and 17. 

[0004] The present invention also provides a method of determining the phase distribution at an image plane, in- 
cluding processing phase-contrast image data of said image. 
15 [0005] The record and therefore the image may be photographic or electronic. The term "image" may thus refer, for 
example, to an observable effect in a set of intensity data, for example a table or other stored record of intensity values; 
the term is not confined to a visual context. The recording medium may comprise a two-dimensional pixilated detector, 
e.g. an electronic detector such as a charge-coupled device (CCD) array. 

[0006] The irradiating means preferably includes a source of x-rays of diameter 20 micron or less, where diameter 
20 refers to the full width of intensity distribution of the source at half maximum intensity. The apparatus may advanta- 
geously further include a suitable stage or holder for samples containing the internal boundary being imaged. 
[0007] The penetrating radiation, e.g. x-ray radiation, may be polychromatic and is preferably in the hard x-ray range, 
i.e. in the range 1 keV to 1 MeV. 

[0008] The separation of the boundary and the detecting means is preferably selected to enhance the resolution of 
25 the image. For example, it has been observed that a sharper image, i.e. one with better contrast, is achieved by in- 
creasing separation. For instance contrast is improved at least for a separation of about 1 m relative to a separation 
of 0.4 m. This may partly be because background noise is diminished with increasing separation but the intensity 
variation effect arising from the change in the local direction of propagation is substantially preserved. 
[0009] The term "lateral spatial coherence" herein refers to the correlation of the complex amplitudes of waves be- 
30 tween different points transverse to the direction of propagation of the waves. Lateral spatial coherence is said to occur 
when each point on a wavefront has a direction of propagation which does not change over time. In practice, high 
lateral spatial coherence may, for example, be achieved by using a source of small effective size or by observing the 
beam at a large distance from the source. For example, for 20 keV x-rays a source size of 20 urn diameter or less 
would typically be appropriate. The smaller the source size the better for the purposes of this invention, provided total 
35 flux from the source is sufficient. Lateral spatial coherence may need to be preserved by careful selection of the x-ray 
window of the source, e.g. such that it is of highly uniform thickness and homogeneity. 

[001 0] Preferred embodiments of the present invention are hereinafter described, by way of example only, with ref- 
erence to the accompanying Figures, in which: 

40 Figure 1 is a diagram, presented in three parts for purposes of illustration, showing a circular cross-section object 

being irradiated by a parallel beam; 

Figure 2 is a diagram of a circular cross-section object being irradiated by a polychromatic beam and the intensity 
of the phase-contrast image produced; 

Figure 3 is a diagram of an x-ray optics configuration according to an embodiment of the invention; and 
45 Figures 4 and 5 are x-ray images of various boundaries derived in accordance with the invention, as subsequently 

detailed herein. 

[0011] It is first now proposed to outline the mathematical basis of the present invention. 

[0012] Variations in thickness and x-ray refractive index, n(X) = 1-8(X)-iP(X), of a sample will invariably lead to a 
50 change in the shape of an x-ray wavefront on passing through the sample. The real component 1-5(A.) of n relates to 
the degree of refraction and the imaginary component -i(3(X) relates to the degree of absorption. More specifically, for 
a single element substance 
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P(X) = A H(X) (2) 

where u.(X) is the linear absorption coefficient, r Q is the classical radius of an electron, N 0 is the number of atoms per 
5 unit volume and f R is the real part of the atomic scattering factor at zero scattering angle. The coefficient 8 is proportional 
to X 2 and (3 is proportional to X 4 and also X is proportional to 1 /energy of the x-ray photon emitted from the source. 
[001 3] The magnitude of the wavefront distortions is related to the gradient of the phase variations transverse to the 
direction of propagation of the wavefront. In the geometrical optics approximation, the phase difference, <j>, for a ray 
path through an object is proportional to the integral of the decrement of the real part of the refractive index, 8, along 
10 that ray path. For the coordinate system illustrated in Figure 1 this can be expressed generally as 
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4>(*,z) = */* {nix**) - 1) dz' 0) 



where k is equal to 2tc/X. The angular deviation Act of the local scattered wavevector from that of the local incident 
wavevector is proportional to the gradient of the phase difference in the direction perpendicular to the local incident 
wavevector. The word "local" refers to a point (x,y,z) on the wavefront. Mathematically the local scattered wavevector 
can be written for the coordinate system illustrated in Figure 1 as 



where s(x,y,z) is the normal to the wavefront at point (x.y.z) and the above relationship is valid in the paraxial approx- 
imation when (3<t>/3x) 2 + (dtyfdy) 2 « k 2 . The angular deviation Act can be expressed as 

A « - ±**™ = n {a&g . il A / (s) 

k ax j « I 3r J 



[0014] The angular deviation Act is therefore dependent on a refractive index variation perpendicular to a propagation 
wavevector, and the amount of deviation depends on the length over which the variation occurs in the direction of the 
wavevector, e.g. the thickness of a sample. 
40 [001 5] To illustrate the nature of the effect, consider the case of a spherical object, n, of refractive index n M embedded 
in a medium of refractive index n Q = 1 , as illustrated in Figures 1 and 2. 

[0016] The x-ray optical path length differences through the sample relative to through a vacuum lead to a phase 
difference <J>(x) and hence to a phase gradient d<j>/3x in the direction (Figure 1 ) transverse to the local direction of 
propagation. The phase difference between ray 1 which passes through the object parallel to the z-axis at constant 
45 distance from it and the reference ray 0 is: 



4>(x.y) = HM - 22- 6<\)z(xj). (6) 

where z(x,y) is the length of the intersection of ray 1 with CI and 

z(x,y) = 2*jR 2 -x 2 -y 2 , (7) 
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and R is the radius of Q and 5 is the decrement of refractive index coefficient. Mathematically, for a circular sectioned 
object in the x-z plane, the expression for d$/dx and the angular deviation Aa between an incident ray and the corre- 
sponding refractive ray for a given x is: 



10 



4 

K 

2n ax 



(8) 
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In equation (8), 5(A.) is slowly varying and it can readily be seen that the phase gradient diverges at x = ±R, where the 
rays can deviate by very large angles from the optic axis. In these limits, the angular deviations of the scattered beams 
can be very large and lead to an observable loss in intensity I in the corresponding forward direction, the position of 
which is independent of wavelength, as shown in Figure 2 for a polychromatic beam B. The decrement of refractive 
index coefficient, S, is typically of order 10" 5 to 10" 6 for light elements and hard x-rays but nevertheless the deviation 
angle Aa may be quite large when x is close to ±R, i.e. at the boundary of the sample or at an internal boundary feature. 
[0017] The nature of the contrast obtained under different conditions of source size, object-source distance and 
object-image distance, and also the spectral distribution of the source need to be considered. A further consideration 
affecting contrast is the degree of modification of the wavefront introduced by the object. 

[001 8] For the plane-wave case, to help understand the role of these factors on contrast in image formation, we can 
to a first approximation use the formula derived by Cowley (J.M. Cowley, "Diffraction Physics", 2nd Ed., p.60, North 
Holland, 1981 ) for the Fresnel diffraction contrast from a phase object. According to this formal, for a one-dimensional 
phase object producing a phase change, <t>(x), under plane-wave illumination with wavelength X, the intensity distribution 
at a distance R 2 from the object is given by 



R 2 X 

/(*>=1 + -£-♦"<*) 



(9) 



35 which is valid to first order in the quantity (R 2 XI2n) $"(*). assumed small. From this apparently simple formula, one 
can draw some significant conclusions, namely: 

i) the contrast varies directly with R 2 , 

ii) the structure of the image is ^-independent. Only the contrast is affected. For a polychromatic source one would 
40 simply replace X in the formula by a spectrally weighted sum. 

[0019] To get some feeling for the range of validity of the above formula for the present x-ray case, let us suppose 
there is an object feature for which the phase transmitted by the object varies by 1 radian over a lateral distance of 1 0 
microns. Then <J>" — 10 10 nr 2 , and for X ~ 1A, R 2 ~ 1m, we see that (R 2 X/2k) '"(x) < 1. Thus the formula should be 

45 valid even for small phase objects or reasonably rapid variations in phase. However, for very sharp edges or changes 
of slope, such as are often used in calculations of artificial test objects (e.g. fibres), may become too large (even 
infinite), so the formula breaks down. But even in these cases the general form of the image (a black/white line from 
a sharp step object) is reproduced but not the subsidiary fringes typical of diffraction from such discontinuities. On the 
other hand, and probably of more practical significance, we see that for smaller <j>"(x), i.e. larger features with less rapid 

50 lateral variation, the contrast will be low, and may well limit the practical visibility. 

[0020] A more exact mathematical treatment of this type of imaging with plane-waves has recently been carried out 
in terms of Fresnel diffraction (P. Cloetens, R. Barrett, J. Baruchel, J. P. Guigay and M. Schlenker, J. Phys. D.: Appld. 
Phys., 1996 29, 133-46; J. P. Guigay, Optik, 1977 49, 121-5). Their treatment gives the same equation as presented 
above to first order. However the more accurate treatment leads to the conclusion that the maximum contrast for a 

55 spatial frequency u occurs when 2 X R 2 u 2 = 1 , at least for the normal range of conditions expected in phase-contrast 
radiography. The spatial frequency u relates to the structure of the object being imaged, where u equals 1/A where A 
is the spatial period of a Fourier component of the imaged object. 

[0021] These treatments all refer to illumination with an ideal plane-wave. Any divergence in the beam will blur the 
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image by an amount proportional to R 2 (in this respect behaving in the same way as in conventional radiography). The 
above authors (Cloetens et al.) then show that the overall optimum R 2 , taking into account both contrast and resolution 
is given by 

R 2 < 2 XI a (10) 



where a = sIR^ is the angle subtended by the source at the object and relates to the (almost) plane-wave case. It 
should be noted that Cloetens et al. specifically prescribe the need for a highly monochromatic source of x-rays and 
consider only the plane-wave case in contradistinction to the preferred embodiments described herein. 
[0022] As pointed out, the treatments above relate specifically to the plane-wave case whereas we are principally 
concerned with the spherical-wave case which more closely relates to convention radiography. To help understand the 
spatial-wave case, we now consider the relationship between the two which can be usefully established via a simple 
analysis of the Fresnel-Kirchhoff expression for imaging an object with a point source at a distance R-j from the object 
(the spherical-wave case). This shows that there is a simple relationship for the spherical-wave case involving terms 
of the plane-wave case but with a modified object to image distance R', such that 



l=± + ± (11) 

and with the image magnified by (R^ + R 2 )/R v From simple geometrical arguments, based on ray optics, it appears 
that loss of contrast or resolution due to source size will not be a problem in the spherical-wave case as both the image 
and the source size are magnified, the latter by R 2 /Ri which asymptote to the same factor for large R 2 . The factor 

25 affecting contrast for the spherical-wave case is that (for the range of energies and spatial resolutions relevant for 
radiography) 2 X R 2 (1 + R 2 /R<|) u 2 should be large (but is typically less than 1 ). This expression may be large because 
R 2 is large, or X large or the spatial frequency, u, is large. As an illustration, for practical radiographic purposes the 
following values might serve as being indicative X = 0.2 A; u < 2 x 10 5 (corresponding to a spatial period of 20 micron 
or more) so that R 2 ~ 2.5 m (assuming R 2 IR^ = 3, say) would give maximum contrast for the highest spatial frequency. 

30 Larger values of R 2 would be appropriate for maximum contrast of lower spatial frequencies. 

[0023] It may be noted that the function will tend to enhance the edges and boundaries of a phase object in an 
image. If there is also an absorptive component of the object it will, at least to first order, add directly to the image 
contrast (e.g. see Equation 7 in Guigay, 1977). The present technique could complement and enhance the usual 
radiological image, as well as yielding new information. It may also be noted that proper treatment of the contrast in 

55 an image involving differential phase-contrast (involving the Laplacian of <(>) requires numerical processing of the image 
via, for example, solution of the transport of intensity equations (see T.E. Gureyev, A. Roberts and K. Nugent, Journal 
of Optical Society of America, Vol. A12, pp.1932 and pp.1942, 1995 incorporated herein by reference) in order to 
retrieve the phase, <|>(x). 

[0024] We turn now to practical arrangements for applying the concept arising from these determinations. In a first 
40 embodiment (Figure 3), there is a source S of high spatial coherence and an x-ray imaging detector D, for example 
film, photostimulable phosphor plates (e.g. Fuji Image Plates), or a two-dimensional electronic detector. Regions of 
sharp refractive index variation transverse to the direction of propagation, or thickness variation in the direction of 
propagation, can lead to a significant change in the local direction of propagation of the wavefront passing through 
those regions. Thus a spherical wavefront W1 emanating from the point source S becomes distorted to W2 on passing 
45 through the object O. By recording the intensity of the wavefront at a sufficient distance from the sample, intensity 
variations due to sharp refractive index and thickness variations in the sample may be detected and their location 
recorded in an image. This corresponds to a form of differential phase-contrast imaging. The location of the imaging 
detector is chosen such that the spatial resolution of the detector is sufficient to resolve the intensity differences arising 
from the severe distortions of the wavefront and to optimise contrast, as described above, subject to practical consid- 
50 erations. 

[0025] Typically, the sharp gradients in refractive index or thickness will be imaged as sharp losses or rapid variation 
in intensity at corresponding points in the image. This feature of intensity loss or rapid variation at a given point in the 
image is essentially independent of wavelength and can therefore lead to very sharp contrast variations in the image 
even when a polychromatic source is used. 
55 [0026] This configuration has the feature that for a circular source distribution, the spatial resolution in the image is 
the same for both directions and is essentially determined by the source size. It also has the advantage that considerable 
magnification of the image is possible and so recording media such as Fuji Image Plates may be used which have 
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many desirable properties such as wide dynamic range and high sensitivity but not high spatial resolution. 
[0027] In addition to the source and detector involved in this configuration, a high resolution angular analyser may 
be inserted between the sample and the detector. The high resolution angular analyser might for example be a suitably 
curved crystal in Laue geometry with curvature chosen for some appropriate characteristic wavelength of the source. 
5 This variation in the method is aimed at resolving weaker variations in refractive index and thickness of the sample 
than are observable with the first described configuration. 

[0028] It may be noted that a very substantial magnification of the image is possible so that very high spatial resolution 
in the image may be achieved even with much lower spatial resolution detectors such as Fuji Image Plates. Also it 
may be noted that since the method of image formation is essentially independent of x-ray energy, the sources can be 
10 operated at high tube voltage and so lead to lower absorbed dose to the sample, which is important in clinical appli- 
cations. 

[0029] Some examples of phase-contrast images recorded using the aforementioned technique are illustrated in 
Figures 4 and 5. Figure 4 shows an image of the edge of a 10um plastic film which is the same as that used in Davis, 
Gao, Gureyev, Stevenson and Wilkins (Phys. Rev. Letters, 1995, Vol. 74, p. 3173) and corresponds to a pure phase 

15 object. Figure 5 shows images of an air bubble and glass fibres in a polymer matrix based on a similar sample to that 
reported in Davis, Gao, Gureyev, Stevenson and Wilkins (Nature Vol. 373 pp. 595-8, 1995) and corresponds to an 
almost pure phase object. In each case clear additional contrast can be seen over that expected for a normal absorption- 
contrast image. In particular, in Figure 4 the edge of the film is clearly visible as a black/white contrast feature as also 
are the edges of the bubbles and the fibres. The source used was a nominal 10pm diameter microfocus source (Kevex 

20 Model PXS) with Cu anode operated at 40 kV. For Figure 4 the source to sample and sample to film distances were 
both 700 mm while for Figure 5 the corresponding distances were 120 mm and 1000 mm, respectively. It should be 
noted that contrast in the present instances is visible almost entirely due to the high spatial coherence of the source. 
The contrast is primarily an intensity loss contrast and in that sense resembles normal absorption but is different in 
that it represents an intensity loss due to refractive scattering (or Fresnel diffraction) at the object boundaries as shown 

25 by equation (8). A normal fine focus source of diameter 0.1 mm would have a projected size of approximately the length ) 
of the 0.1 mm scale bar shown on the photographs and so largely smear out this contrast. \ 
[0030] To provide a comparison of phase-contrast imaging as described herein and standard absorption imaging, 't 
the table below sets out the absorption thickness tg of a carbon sample required to achieve 65% absorption and the K 
phase thickness tp of the sample required to achieve a phase change in <j> of 2k, for different source energies E. 

30 



TABLE 1 



35 



E keV 


MA) 


t a (nm) 


tp(^m) 


50 


0.25 


435000 


133 


12 


1 


5000 


30 


1.2 


10 


4 


3 


0.25 


50 


1.3 


1.2 



[0031] The results in the table illustrate how phase-contrast imaging can be used to image very small objects with 
high energy sources. 

[0032] Advantageously, the beam path between sample and detector may involve evacuated tubes with x-ray trans- 
parent windows or similar means to reduce the effects of air scattering making sure that their optical quality is such 
that they do not have a detrimental effect on the coherence of the x-ray beam. 

[0033] The present method should be especially well suited to imaging of such features as cracks, voids and dela- 
mination effects in various types of materials, since these features involve maximum differences in x-ray refractive 
index and the spatial variation can be extremely sharp. To give observable contrast, the source is preferably of a very 
small effective size, say less than of order 20u.m, and the detector is preferably a high resolution imaging detector such 
as x-ray film or a two-dimensional electronic detector, e.g. a CCD array. The method may also prove useful in signifi- 
cantly enhancing the contrast of important features in clinical radiography. 

[0034] The present application outlines some simplified conditions and configurations for differential phase-contrast 
imaging using penetrating radiation such as hard x-rays, which are particularly aimed at clinical and industrial appli- 
cations. These new approaches are more closely related to traditional methods used for absorption-contrast radiog- 
raphy and should be easier to implement than our earlier described methods of the aforementioned WO95/05725 and 
PN5811/95, especially for large areas of irradiation. They should also have considerably shorter exposure times for a 
given source power than the earlier monochromatic methods since they may use a very wide spectrum from the source. 
[0035] Throughout this specification , unless the context requires otherwise, the word "comprise", or variations such 
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as "comprises" or "comprising", will be understood to imply the inclusion of a stated integer or group of integers but 
not the exclusion of any other integer or group of integers. 

5 Claims 

1 . A method of deriving a phase-contrast image of an internal boundary of an object (O), which boundary represents 
a refractive index variation, comprising: 

10 irradiating the boundary with a propagated wavefront (W1 ) of x-ray radiation that penetrates the object and 

that has a propagation direction such there is a significant component of the propagation vector transverse to 
the direction of said refractive index variation, and further having a lateral spatial coherence sufficiently high 
for the variation in refractive index to cause a detectable change in the local direction of propagation of the x- 
ray radiation wavefront at the boundary that is manifested as contrast via Fresnel diffraction; and 

15 detecting and recording intensity of at least a portion of said wavefront (W2) of radiation, after it has traversed 

said boundary and emerged from said object, in a manner whereby an effect of said change in the local direction 
of propagation is observable in the intensity record as said contrast and thereby recorded as a local diminution 
or rapid variation of intensity of the x-ray radiation, whereby the intensity record substantially forms an image 
of the boundary. 

20 

2. A method as claimed in claim 1 , wherein said radiation is polychromatic x-ray radiation. 

3. A method as claimed in claim 1 or 2, wherein said radiation is x-ray radiation having energy in the range 1 keV to 
1 MeV. 

25 

4. A method as claimed in any one of claims 1 to 3, including irradiating said boundary with a source of x-rays having 
a diameter transverse to the optic axis less than or equal to 20 urn. 

5. A method as claimed in any one of claims 1 to 4, including separating said boundary and the position of detecting 
30 said portion of said radiation by a distance which enhances the contrast and/or resolution of the part of an image 

comprising the record of said local diminution or rapid variation of wavefront intensity. 

6. A method as claimed in claim 5, wherein said distance is greater than or equal to 0.3 m. 

35 7. A method as claimed in claim 6, wherein said distance is greater than or equal to 0.7 m. 

8. A method as claimed in any one of claims 1 to 7, wherein said variation in the detected intensity of said wavefront 
is sharp and localized. 

40 9. A method of determining the phase distribution at an image plane, including processing a record derived by a 
method as claimed in any one of claims 1 to 8. 

10. A method of determining the phase distribution at an image plane as claimed in claim 9, wherein said processing 
is based on Maxwell's equations for electromagnetic radiation. 

45 

11. A method of determining the phase distribution at an image plane as claimed in claim 10, wherein said equations 
are the transport of intensity equations. 

12. A method as claimed in any one of claims 1 to 9, wherein said radiation is generated by and provided from a source 
50 to said boundary without Bragg diffraction. 

13. A method as claimed in any one of claims 1 to 9 further including processing said record so as to derive from the 
record said local diminution or rapid variation of the detected intensity of said x-ray radiation and so identify the 
representation of the boundary. 



55 



14. A method as claimed in claim 1 3 wherein said step of irradiating said boundary comprises irradiating said boundary 
with an unfocused propagated wavefront of x-ray radiation, and wherein said step of detecting intensity comprises 
detecting intensity of at least a portion of said wavefront of said x-ray radiation passing through said boundary so 
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as to form said record without focusing said wavefront after it passes through said boundary. 

15. A method as claimed in claim 13 or 14, wherein said processing is based on Maxwell's equations for electromag- 
netic radiation. 

16. A method as claimed in claim 15, wherein said equations are the transport of intensity equations. 

17. An apparatus for deriving a phase-contrast image of an internal boundary of an object (O), which boundary rep- 
resents a refractive index variation, comprising: 



means (S) to irradiate the boundary with a propagated wavefront (W1 ) of x-ray radiation that penetrates the 
object and that has a propagation direction such that there is a significant component of the propagation vector 
transverse to the direction of said refractive index variation, and further having a lateral spatial coherence 
sufficiently high for the variation in refractive index to cause a detectable change in the local direction of prop- 
15 agation of the x-ray radiation wavefront at the boundary that is manifested as contrast via Fresnel diffraction; 

and 

means (D) to detect and record intensity of at least a portion of said wavefront (W2) of radiation after it has 
traversed said boundary and emerged from said object, in a manner whereby an effect of said change in the 
local direction of propagation is observable in the intensity record as said contrast and thereby recorded as a 
20 local diminution or rapid variation of intensity of the wavefront of radiation, whereby the intensity record sub- 

stantially forms an image of the boundary. 

18. An apparatus as claimed in claim 17, wherein said radiation is polychromatic x-ray radiation. 

25 19. An apparatus as claimed in claim 17 or 18, wherein said radiation is x-ray radiation having energy in the range 1 
keV to 1 MeV. 

20. An apparatus as claimed in claim 17, 18 or 19, wherein said means to irradiate is a source of x-rays having a 
diameter transverse to the optic axis less than or equal to 20 urn. 

30 

21 . An apparatus as claimed in any one of claims 1 7 to 20, further including holder means to hold an object containing 
said boundary and so locate the boundary at a predetermined position, whereby the separation of said boundary 
and the position of detecting said portion of said radiation may be set at a distance which enhances the contrast 
and/or resolution for part of an image comprising the record of said local diminution or rapid wavefront variation 

35 of intensity. 

22. An apparatus as claimed in claim 21, wherein said detection means and said holder means are disposed so that 
said distance is greater than or equal to 0.3 m. 

40 23. An apparatus as claimed in claim 21, wherein said detection means and said holder means are disposed so that 
said distance is greater than or equal to 0.7 m. 

24. An apparatus as claimed in any one of claims 17 to 23, wherein said radiation is generated by and provided from 
a source to said boundary without Bragg diffraction. 



25. An apparatus according to any one of claims 17 to 24 further including means for processing said record so as to 
derive from the record said corresponding variation in the detected intensity of said wavefront in said record and 
so identify the representation of the boundary. 

Patentanspriiche 

1. Verfahren zum Gewinnen eines Phasenkontrastbilds einer internen Grenze eines Objekts (O), wobei diese Grenze 
eine Brechungsindexabweichung darstellt, das folgendes umfasst: 

Bestrahlen der Grenze mit einer sich ausbreitenden Wellenfront (W1 ) von Rontgenstrahlung, die in das Objekt 
eindringt und die eine Ausbreitungsrichtung hat, so dass eine wesentliche Komponente des Ausbreitungsvek- 
tors quer zur Richtung der Brechungsindexabweichung vorhanden ist und die weiter eine laterale raumliche 
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Koharenz aufweist, die ausreichend hoch ist, dass die Abweichung des Brechungsindex eine erfassbare An- 
derung der ortlichen Richtung der Ausbreitung der Rontgenstrahlungswellenfront an der Grenze verursacht, 
die sich uber Fresnelsche Beugung als Kontrast offenbart; und 

Erfassen und Aufzeichnen von Intensitat von mindestens einem Teil der Wellenfront (W2) der Strahlung, nach- 
5 dem Sie die Grenze uberquert hat und aus dem Objekt ausgetreten ist, derart, dass ein Effekt der Anderung 

der ortlichen Richtung der Ausbreitung in der Intensitatsaufzeichnung als dieser Kontrast beobachtbar ist und 
dadurch als ortliche Abnahme oder schnelle Abweichung der Intensity der Rontgenstrahlung aufgezeichnet 
wird, wodurch die Intensitatsaufzeichnung im Wesentlichen ein Bild der Grenze bildet. 

10 2. Verfahren nach Anspruch 1 , wobei es sich bei der Strahlung um polychromatische Rontgenstrahlung handelt. 

3. Verfahren nach Anspruch 1 oder 2, wobei es sich bei der Strahlung um Rdntgenstrahlung mit Energie im Bereich 
von 1 keV bis 1 MeV handelt. 

15 4. Verfahren nach einem der Anspruche 1 bis 3, das das Bestrahlen der Grenze mit einer Quelle von Rontgenstrahlen 
mit einem Durchmesser quer zur optischen Achse von 20 jxm oder weniger einschliefit. 

5. Verfahren nach einem der Anspruche 1 bis 4, das das Trennen der Grenze und der Position der Erfassung der 
Strahlung um einen Abstand einschlie&t, der den Kontrast und/oder die Auflbsung des Teils eines Bilds verbessert, 

20 der die Aufzeichnung der ortlichen Abnahme beziehungsweise schnellen Abweichung der Wellenfrontintensitat 

umfasst. 

6. Verfahren nach Anspruch 5, wobei der Abstand grofier oder gleich 0,3 m ist. 

25 7. Verfahren nach Anspruch 6, wobei der Abstand grofier oder gleich 0,7 m ist. 

8. Verfahren nach einem der Anspruche 1 bis 7, wobei die Abweichung in der erfassten Intensitat der Wellenfront 
scharf und ortlich begrenzt ist. 

30 9. Verfahren zur Bestimmung der Phasenverteilung an einer Bildebene, das die Verarbeitung einer durch ein Ver- 
fahren nach einem der Anspruche 1 bis 8 gewonnenen Aufzeichnung einschlie&t. 

1 0. Verfahren zur Bestimmung der Phasenverteilung an einer Bildebene nach Anspruch 9, wobei die Verarbeitung auf 
den Maxwellschen Gleichungen fur elektromagnetische Strahlung beruht. 

35 

11. Verfahren zur Bestimmung der Phasenverteilung an einer Bildebene nach Anspruch 10, wobei es sich bei den 
Gleichungen um die Intensitatstransportgleichungen handelt. 

12. Verfahren nach einem der Anspruche 1 bis 9, wobei die Strahlung ohne Braggsche Beugung von einer Quelle 
*o erzeugt und von dort an die Grenze geliefert wird. 

13. Verfahren nach einem der Anspruche 1 bis 9, das weiter das Verarbeiten der Aufzeichnung einschlielit, so dass 
die ortliche Abnahme beziehungsweise die schnelle Abweichung der erfassten Intensitat der Rontgenstrahlung 
aus der Aufzeichnung gewonnen und so die Darstellung der Grenze identifiziert wird. 

45 

14. Verfahren nach Anspruch 13, wobei der Schritt des Bestrahlens der Grenze das Bestrahlen der Grenze mit einer 
nicht fokussierten sich ausbreitenden Wellenfront von Rontgenstrahlen umfasst und wobei der Schritt des Erfas- 
sens von Intensitat das Erfassen von Intensitat von mindestens einem Teil der durch die Grenze laufenden Wel- 
lenfront der Rontgenstrahlung umfasst, um die Aufzeichnung zu bilden, ohne die Wellenfront zu fokussieren, nach- 

50 dem sie die Grenze durchlaufen hat. 

15. Verfahren nach Anspruch 13 oder 14, wobei die Verarbeitung auf den Maxwellschen Gleichungen fur elektroma- 
gnetische Strahlung beruht. 

55 16. Verfahren nach Anspruch 15, wobei es sich bei den Gleichungen um die Intensitatstransportgleichungen handelt. 

17. Vorrichtung zum Gewinnen eines Phasenkontrastbilds einer internen Grenze eines Objekts (O), wobei diese Gren- 
ze eine Brechungsindexabweichung darstellt, die folgendes umfasst: 
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ein Mittel (S) zum Bestrahlen der Grenze mit einer sich ausbreitenden Weilenfront (W1 ) von Rontgenstrahlen, 
die in das Objekt eindringt und die eine Ausbreitungsrichtung hat, so dass eine wesentliche Komponente des 
Ausbreitungsvektors quer zur Richtung der Brechungsindexabweichung vorhanden ist und die weiter eine 
laterale raumliche Koharenz aufweist, die ausreichend hoch ist, dass die Abweichung des Brechungsindex 
5 eine erfassbare Anderung der ortlichen Richtung der Ausbreitung der RSntgenstrahlungswellenfront an der 

Grenze verursacht, die sich uber Fresnelsche Beugung als Kontrast offenbart; und 

ein Mittel (D) zum Erfassen und Aufzeichnen von mindestens einem Teil der Weilenfront (W2) der Strahlung, 
nachdem sie die Grenze uberquert hat und aus dem Objekt ausgetreten ist, derart, dass ein Effekt der Ande- 
rung der ortlichen Richtung der Ausbreitung in der Intensitatsaufzeichnung als dieser Kontrast beobachtbar 
10 ist und dadurch als ortliche Abnahme Oder schnelle Abweichung der Intensitat der Rontgenstrahlung aufge- 

zeichnet wird, wodurch die Intensitatsaufzeichnung im Wesentlichen ein Bild der Grenze bildet 

18. Vorrichtung nach Anspruch 17, wobei es sich bei der Strahlung um polychromatische Rontgenstrahlung handelt. 

15 19. Vorrichtung nach Anspruch 17 oder 18, wobei es sich bei der Strahlung um Rontgenstrahlung mit Energie im 
Bereich von 1 keV bis 1 MeV handelt. 

20. Vorrichtung nach Anspruch 17, 18 oder 19, wobei es sich beim Mittel zum Bestrahlen um eine Quelle von Ront- 
genstrahlen mit einem Durchmesser quer zur optischen Achse von 20 \im oder weniger handelt. 

20 

21. Vorrichtung nach einem der Anspriiche 17 bis 20, die weiter ein Haltemittel zum Halten eines die Grenze enthal- 
tenden Objekts einschliefit, um so die Grenze an einer vorbestimmten Position zu positionieren, wodurch fur die 
Trennung der Grenze und der Position der Erfassung des Teils der Strahlung ein Abstand eingestellt werden kann, 
der den Kontrast und/oder die Auflbsung fur etnen Teil eines Bilds verbessert, der die Aufzeichnung der ortlichen 

25 Abnahme beziehungsweise schnellen Abweichung der Intensitat umfasst. 

22. Vorrichtung nach Anspruch 21, wobei das Erfassungsmittel und das Haltemittel so angeordnet sind, dass der 
Abstand grower oder gleich 0,3 m ist. 

30 23. Vorrichtung nach Anspruch 21, wobei das Erfassungsmittel und das Haltemittel so angeordnet sind, dass der 
Abstand grofter oder gleich 0,7 m ist. 

24. Vorrichtung nach einem der Anspruche 17 bis 23, wobei die Strahlung ohne Braggsche Beugung von einer Quelle 
erzeugt und von dort an die Grenze geliefert wird. 

35 

25. Vorrichtung nach einem der Anspruche 17 bis 24, die weiter ein Mittel zum Verarbeiten der Aufzeichnung ein- 
schliefM, um die entsprechende Abweichung in der erfassten Intensitat der Weilenfront aus der Aufzeichnung zu 
gewinnen und so die Darstellung der Grenze zu identifizieren. 

40 

Revendications 

1. Precede de derivation d'une image a contraste de phase d'une limite interne d'un objet (O), limite qui represente 
une variation d'indice de refraction, qui comprend : 

45 

I'irradiation de la limite avec un front d'onde propage (W1) d'un rayonnement X qui penetre dans I'objet et a 
un sens de propagation tel qu'il y a un element significatif du vecteur de propagation transversal par rapport 
au sens de ladite variation d'indice de refraction, et qui a en outre une coherence spatiale laterale sufflsamment 
elevee pour que la variation de I'indice de refraction provoque un changement detectable dans le sens local 
50 de la propagation du front d'onde de rayonnement X a la limite qui se manifeste en tant que contraste par la 

diffraction de Fresnel ; et 

la detection et I'enregistrement de I'intensite d'au moins une portion dudit front d'onde (W2) de rayonnement, 
apres que celui-ci ait traverse ladite limite et soit sorti dudit objet, d'une maniere par laquelle un effet dudit 
changement du sens local de propagation est observable dans I'enregistrement d'intensite comme etant ledit 
55 contraste et par consequent est enregistre comme diminution locale ou variation rapide de I'intensite du rayon- 

nement X, par lequel I'enregistrement d'intensite forme en grande partie une image de la limite. 

2. Procede tel que revendique a la revendication 1 , dans lequel ledit rayonnement est un rayonnement X polychro- 
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matique. 

3. Procede tel que revendique a la revendication 1 ou 2, dans lequel ledit rayonnement est un rayonnement X pos- 
sedant une energie de I'ordre de 1 keV a 1 MeV. 

5 

4. Precede tel que revendique dans Tune quelconque des revendications 1 a 3, qui incorpore I'irradiation de ladite 
limite avec une source de rayons X ayant un diametre transversal par rapport a I'axe optique inferieur ou egal a 
20 \im. 

10 5. Procede tel que revendique dans Tune quelconque des revendications 1 a 4, qui incorpore la separation de ladite 
limite et la position de detection de ladite portion dudit rayonnement par une distance qui ameliore le contraste et/ 
ou la resolution de la partie d'une image qui comprend I'enregistrement de ladite diminution locale ou variation 
rapide de I'intensite de front d'onde. 

15 6. Procede tel que revendique a la revendication 5, dans lequel ladite distance est superieure ou egale a 0,3 m. 

7. Procede tel que revendique a la revendication 6, dans lequel ladite distance est superieure ou egale a 0, 7 m. 

8. Procede tel que revendique dans Tune quelconque des revendications 1 a 7, dans lequel ladite variation dans 
20 I'intensite detectee dudit front d'onde est nette et localisee. 

9. Procede de determination de la distribution de phase a un plan d'image, qui incorpore le traitement d'un enregis- 
trement derive par un procede tel que revendique dans Tune quelconque des revendications 1 a 8. 

25 10. Procede de determination de la distribution de phase a un plan d'image tel que revendique a la revendication 9, 
dans lequel ledit traitement est base sur les equations de Maxwell pour le rayonnement electromagnetique. 

11. Procede de determination de la distribution de phase a un plan d'image tel que revendiqu6 a la revendication 10, 
dans lequel lesdites equations sont le transport d'equations d'intensite. 

30 

12. Procede tel que revendique dans Tune quelconque des revendications 1 a 9, dans lequel ledit rayonnement est 
genere par une source et prevu depuis une source vers ladite limite sans la diffraction de Bragg. 

13. Procede tel que revendique dans Tune quelconque des revendications 1 a 9 qui incorpore en outre le traitement 
35 dudit enregistrement de facon a deriver depuis I'enregistrement ladite diminution locale ou variation rapide de 

I'intensite detectee dudit rayonnement X et a identifier ainsi la representation de la limite. 

14. Procede tel que revendique a la revendication 13 dans lequel ladite etape d'irradiation de ladite limite comprend 
I'irradiation de ladite limite avec un front d'onde propage non focalise d'un rayonnement X, et dans lequel ladite 

40 etape de detection de I'intensite comprend la detection de I'intensite d'au moins une portion dudit front d'onde 

dudit rayonnement X traversant ladite limite de facon a former ledit enregistrement sans focal iser ledit front d'onde 
apres que celui-ci ait traverse ladite limite. 

15. Procede tel que revendique a la revendication 13 ou 14, dans lequel ledit traitement est base sur les equations 
45 de Maxwell pour un rayonnement electromagnetique. 

16. Procede tel que revendique a la revendication 15, dans lequel lesdites equations sont le transport des equations 
d'intensite. 

50 17. Appareil de derivation d'une image a contraste de phase d'une limite interne d'un objet (O), limite qui represente 
une variation d'indice de refraction, qui comprend : 

un moyen d'irradiation de la limite avec un front d'onde propage (W1) d'un rayonnement X qui penetre dans 
I'objet et a un sens de propagation tel qu'il y a un element significatif du vecteur de propagation transversal 
55 par rapport au sens de ladite variation d'indice refraction, et qui a en outre une coherence spatiale laterale 

suffisamment elevee pour que la variation de I'indice refraction provoque un changement detectable dans le 
sens local de propagation du front d'onde de rayonnement X a la limite qui se manifeste en tant que contraste 
par la diffraction de Fresnel ; et 
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un moyen (D) de detection et d'enregistrement de I'intensite d'au moins une portion dudit front d'onde (W2) 
de rayonnement apres qu'il ait traverse ladite limite et soit sorti dudit objet, d'une maniere par laquelle un effet 
dudit changement dans le sens local de propagation est observable dans I'enregistrement d'intensite comme 
etant ledit contraste et par consequent est enregistre comme une diminution locale ou une variation rapide 
d'intensite du front d'onde de rayonnement, par lequel I'enregistrement d'intensite forme en grande partie une 
image de la limite. 

18. Appareil tel que revendique a la revendication 17, dans lequel ledit rayonnement est un rayonnement X polychro- 
matique. 

19. Appareil tel que revendique a la revendication 17 ou 18, dans lequel ledit rayonnement est un rayonnement X 
possedant une energie de I'ordre de 1 keV a 1 MeV. 

20. Appareil tel que revendique a la revendication 17, 18 ou 19, dans lequel ledit moyen d'irradiation est une source 
de rayons X ayant un diametre transversal par rapport a I'axe optique inferieur ou egal a 20 um. 

21. Appareil tel que revendique dans Tune quelconque des revendications 17 a 20, qui incorpore en outre un moyen 
de support d'un objet contenant ladite limite et localiser ainsi la limite a une position predeterminee, par lequel la 
separation de ladite limite et la position de detection de ladite portion dudit rayonnement peuvent etre definies a 
une distance qui ameliore le contraste et/ou la resolution pour une partie de I'image qui comprend I'enregistrement 
de ladite diminution locale ou variation de front d'onde rapide de I'intensite. 

22. Appareil tel que revendique a la revendication 21 , dans lequel ledit moyen de detection et ledit moyen de support 
sont disposes de facon a ce que ladite distance soit superieure ou egale a 0,3 m. 

23. Appareil tel que revendique a la revendication 21, dans lequel ledit moyen de detection et ledit moyen de support 
sont disposes de facon a ce que ladite distance soit superieure ou egale a 0,7 m. 

24. Appareil tel que revendique dans Tune quelconque des revendications 17 a 23, dans lequel ledit rayonnement est 
genere par une source et prevu depuis une source vers ladite limite sans diffraction de Bragg. 

25. Appareil selon Tune quelconque des revendications 1 7 a 24 qui incorpore en outre un moyen de traitement dudit 
enregistrement de facon a deriver depuis I'enregistrement ladite variation correspondante dans I'intensite detectee 
dudit front d'onde dans ledit enregistrement et a identifier ainsi la representation de la limite. 
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